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Abstract 

The accuracy and resource consumption of the four different turbulence models 
based on the eddy viscosity assumption, namely, k — e, two k — cu and Spallart- 
Allmaram models, in modeling airships are investigated. The test airship shape 
is a conventional shape. Three different angles of attack are considered. The 
results are checked against the wind tunnel experimental data. The resource 
consumption study is based on the benchmark of 1500 iterations. Based on all 
data obtained it is evident that Spallart-Allmaras model is the most optimal 
one in the majority of cases. 
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1. Introduction 

Computational Fluid Dynamics allows scientists and engineers to test new air- 
ship shapes for their aerodynamic behaviour using only a computer. Airship 
simulations usually mean employing the flows with high Reynolds numbers. 
This fact requires using turbulence models to predict airship behaviour accu- 
rately. 

The first interest in turbulence modeling appeared before the first computer 
was created, namely, in the beginning of last century. One of the pioneers was 
Prandtl with his mixing-length hypothesis published in 1925 [lj. It was far from 
the modern models, but because all calculations were carried out by hands, the 
main focus was to reduce the number of operations as many as possible. 

The first computers became available for the purposes of scientific research right 
after the Second World War. A new interest in turbulence modeling appeared 
at the same period of time because of the development of jet engines, supersonic 
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aircrafts and some other technologies which required more accurate simulations. 
Many different turbulence models were developed during the period of 1940s- 
1960s. These were the first attempts of accurate prediction of near-wall layer 
turbulence flows. 

But it was the beginning of 1970s when the modern turbulence models started 
to appear. The main achievement was the development of the parent 3 equation 
model by Hanjalic and Launder Q and then the original 2 equation k — e model 
by Launder and Spalding Q. The limitations of the latter model were soon 
found, in particular, inaccurate prediction of low Reynolds near- wall flows. The 
first modification oik — e model for a specific type of flow Q appeared in 1972, 
way before the paper on the finalised original model 0] was published. Other 
turbulence models for the accurate prediction of the boundary layer behaviour 
were developed at the same time (e.g. Ng and Spalding model [f| for turbulent 
kinetic energy k and turbulent length scale I), but k — e with its modification 
became one of the most widespread models in the CFD world. 

In 1980 another iconic turbulence model appeared. It is based on the same 
Boussincsq Hypothesis (or eddy viscosity assumption) and employs the same 
turbulent kinetic energy. But instead of dissipation rate e specific dissipation 
rate to is used. The model was first introduced by Wilcox [6J. Later Menter 
developed a modificatcd model called Menter SST k — uj model 0] , which is 
used along with the original model and lots of its other modifications. 

Another widespread turbulence model called Spallart-Allmaras model was in- 
troduced by Spalart and Allmaras in 1992 Q. This model directly employs one 
equation for turbulent viscosity to close the Reynolds stress tensor in RANS. 
This model was specifically developed for external aerodynamic flows and thus 
is suitable for modeling airships. 

So, in the beginning of 1990s many models and their modifications were available 
to simulate turbulent flows. But the choice of the model for some particular 
applications is not an obvious decision and usually is subject to a separate study. 
That was the main reason for many papers on turbulence models comparison. 
The cases considered in such papers can be divided into two types: developing 
flows and fully developed flows. In the first case strictly transient simulations 
can be used to see the dynamic effects. In the second type the fully developed 
flow can be studied using either transient or steady-state simulations. 

The developing (or unsteady) flows are more complicated and, thus, much more 
papers have been published in this subject. Among the possible applications 
the following studies have been carried out: flows around cube or square (in 2D 
case) jl-Gll 

and square cylinder [12J, simulations of air jets [13J and simulation 



of a working blood pump [14j |. 

The problems related to fully developed flows include, for example, pipe flows 



151 ] and backward facing step problem [16( . The latter is considered as a bench 



mark for turbulence models. One particular interesting application in terms of 
fully developed flows is flow over airships. Omari et.al. [171 ] did the comparison 
between RANS k — e, LES Smagorinsky and VMS-LES models. They used a 
large prolate spheroid shape (airship like shape) with the angle of attack at 20°. 

Modern increasing interest in airships requires better understanding of their 
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aerodynamics and behaviour in order to make them safer and more efficient. 
Such tasks are solved using CFD nowadays. But as it was indicated before, 
turbulence model choice is not a trivial task. Thus a series of special compari- 
son studies for different cases is necessary to give a guidance for scientists and 
engineers in this choice. 

This work is aiming to solve this problem and help with such a choice for the 
case of fully developed flows over the airships for small, medium and large angles 
of attack. In this paper the performance of four turbulence models which work 
with eddy viscosity assumption, namely, two layer realizable k — e, Standard and 
SST Mentor k — u and Spallart-Allmaras turbulence models, are compared. The 
results are validated against the experimental data obtained using wind tunnel. 
Commercial CFD code Star-CCM+ 6.04 was used for all the simulations. 



2. Theory 



All four models used for comparison are based on the Reynolds- Averaged Navier- 
Stokcs equations. Classic form of this equation is 
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where Xi are coordinates, Uj are respective time-averaged velocity components, 
u\ are respective fluctuations of the velocity components, p is a time-averaged 
pressure, p is a density, /j are external forces components, /i is a dynamic 
viscosity. 



The unclosed terms —pu^u'j are called Reynolds stresses. They are caused by 



the fluctuations of the velocity field. Turbulence models are used to close these 
terms. All four models considered in this paper are based on the eddy viscosity 
assumption, which can be written as follows 
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where Sij = ^ {^r- + is a mean rate of strain tensor, v t 

turbulent eddy viscosity, k = and Sij is the Kronecker delta. 
The following four models are used for the study. 



pt/p is a 



2.1. RANS Realizable Two Layer k — e model 

The standard k — e model consists of two transport equation 
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for turbulent kinetic energy k and dissipation rate e respectively. In the equa- 
tions and (01 Pf. - production of k term, Pb - effect of buoyancy term, Y m - 
contribution of dilatation fluctuation to the overall dissipation rate, Sk and S £ 
- user-defined source terms. The constants presented in the equations © and 
((4]) have the following values: 

C u = 1.44, C 2e = 1.92, C M = 0.09, a k = 1.0, o e = 1.3. 

It is not completely clear how to calculate the coefficient C*3 e . The following 
formula has been employed 



C 



3e 



tanh 



\u b \ 



where Vb is the velocity component parallel to the gravitational vector and Ub 
is the velocity component perpendicular to the gravitational vector. 



The turbulent viscosity is represented by the formula 

fit = pGfj, — . 



(5) 



However, the model used in this paper is a combination of realizable modification 
of standard model with two-layer approach. The realizable modification was 
introduced by Shih jlTsj ] . In this modification there is a new transport equation 
for the turbulent dissipation rate e 
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and the coefficient in the model for the turbulent viscosity becomes variable 
and now depends on the mean flow and turbulence properties. 

Two layer approach introduced by Rodi (l9| divides the whole domain into two 
layers. Turbulent dissipation rate s and turbulent viscosity p, are represented as 
a function of the wall distance in the layer adjacent to the wall. The values of 
e are computed using transport equation far from wall and blended smoothly 
with the near-wall function. There is no division into layers for the transport 
equation for the turbulent kinetic energy k, i.e. it is solved in the entire domain. 



2.2. RANS Standard k - uj model 

The model was first introduced by Wilcox Q in 1988. Later he revised this 
model, but there is no validation for complex flows for revised version of the 
model. So, the original model is still widely used and referred to as the standard 
k — uj model. 

This is a two equation model and the main alternative to k — e model. Turbulent 
kinetic energy k and specific dissipation rate to are computed from two transport 
equations and used for turbulent viscosity modeling using the formula 
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The equation for the kinetic energy k is 



dpk dpujk 



dt dxj 
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(8) 



while the equation for the specific dissipation rate uj looks as follows 



dpui dpujU 




dt dxj 



(9) 



where 



e = P*uik 

and constants presented in have the following values 
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This model works very well in viscous boundary layer (even in complicated 
cases) without any modifications. This is a significant advantage for the appli- 
cation of this model in the airship simulations. It also has a disadvantage in 
terms of the frcestream and inlet to value sensitivity. This disadvantage mostly 
affects internal flows, but still it can give some negative effect in case of airship 
simulations as well. 

2.3. RANS SST Menter k - uj model 

The main disadvantage of the standard k — uj model was addressed by Mentor 
in Q, where he introduced his modification for the standard model which was 
supposed to solve the main problem of the model leaving all its advantages. The 
approach is that standard k — uj model is used in the near-wall layer while the 
far field is computed using k — e model with a special function, which blends 
the result of the two models together. There are other minor modifications to 
the standard model in SST model, for example, the last term in the transport 
equation for the specific dissipation rate w becomes 2(1 — Fi){pcr u} 2 / ^)^- ■ 
Again, the details can be found in the original Menter's paper Q. 

2.4- RANS Standard Spallart-Allmaras model 

The Spallart-Allmaras model uses a single transport equation for modified tur- 
bulent viscosity v 



dv du-iV 

i i 
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The turbulent viscosity is modelled using the following formula 

Vt = vf v l, (11) 

where /„i = x 3 / (x 3 + C%i) and x '■— vjv. The notation and constant values for 
equation (|10j) can be found in |Appcndix A| 

This model was specifically developed for the aerospace industry, which gives 
some hopes for accurate and not so resource-intensive (than the other three 
models) solution. Nevertheless, the standard model is usually considered as 
Low Reynolds number flow model, even the authors of the original model have 
shown that it works well for attached flows and flows with mild separation, 
which is exactly the case of airship simulations. That is the main reason for 
choosing the standard model over modified high Reynolds number model. 

3. Experimental and simulation approach 

3.1. Airship shape and environment parameters 

The shape of the model for all the simulations is a classic Zeppelin shape with 4 
fins, all at a 45° to the vertical and horizontal longitudinal planes (see Fig. [TJ. 
The fineness ratio of the airship is 4:1 and the total chord length is 1 m. The 
centre of gravity is located on the axis of symmetry of the airship, the distance 
from the nose is 0.451 m. The simulations were conducted for models with for 
three angles of attack, namely —0.4°, 11.62° and 35.62°. 
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Figure 1: The model shape (dimensions are in mm) 

The ambient medium is a gas with the following parameters: 

• constant density: p = 1.204 kg/m 3 , 

• dynamic viscosity: fj, = 1.789 ■ 10~ 5 Pa ■ s, 

• ambient pressure at the sea level: po = 101325 Pa, 

• airship velocity: v = 37 m/s. 
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The coordinate system originates at the centre of gravity and consists of three 
axes X, Y and Z fixed to the body of the airship. The origin is in the centre 
of gravity. X represents the axis of symmetry of the airship looking backwards, 
Z is an axis perpendicular to X and looking up, Y makes up a right-handed 
coordinate system. 

The air domain was chosen as a block and in cross section had the same sizes as 
the wind tunnel cross section in the test zone. The details about the wind tunnel 
tests see in Section l3~2l The boundary conditions again simulate the wind tunnel 
behaviour but with some modifications. The front surface of the air domain is 
chosen as the velocity inlet and the velocity was set to v = 37 m/s. The back 
surface was simulating the outlet with ambient pressure po — 101325 Pa. The 
sidcwalls of the block were made as slip surfaces and the surface of the airship 
itself was made as non-slip one. 

The results are presented as the aerodynamic coefficients rather than the forces 
and moments. For the coefficients, the following reference values were chosen: 

• reference density: p = 1.204 kg/m 3 , 

• reference area: S = 0.101 to 2 , 

• reference length: I = 0.318to, 

• reference velocity: v = 37 m/s. 

All the simulations arc of the steady-state type. We required each simulations 
to make 1500 iterations. 

3.2. Experimental setup 

The experimental data was obtained using a closed wind tunnel with a 2x2 
meters cross section in the test area. The airship sizes and the flow velocity was 
the same as for the simulation. The airship model was made out of polystyrene 
foam and fixed in supposed centre of gravity (see above). The sizes and the 
shape of the model are exactly the same as can be seen in the Figure \T\ All the 
experiments were duplicated and an average of two values was taken as final 
result. 

3.3. Mesh parameters 

A polyhedral mesh is a usual choice for the aerospace application. In particular, 
airships have quite steep filleted or rounded surfaces and polyhedral mesh is 
the best type of mesh for such shapes. Near-wall prism mesh consisting of 10 
layers was used for more accurate simulation of the boundary layer flows. The 
air domain for each simulation contained approximately 4.5 million cells. The 
mesh is built in such a way that Y+ values in the nearest to the wall surface 
cells are around unity. 
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4. Results and discussion 



The results are presented in this section in Tables [T][3] for small, medium and 
large angles of attack. Each table contains two force coefficients c x and c z (the 
forces in the directions of the axes X and Z respectively) from wind tunnel test 
and CFD simulations using 4 different turbulence models described in section [2] 
together with the errors in comparison with experimental data. CPU time used 
for doing the benchmark of 1500 iterations can also be found in the tables as 
well as on the Figure [5] Comparison between experiment results and modeling 
solutions using all 4 turbulence models arc presented in a graphical form in 
Figures [3] and [U 

cpu time, s 




Figure 2: CPU time required for doing a benchmark of 1500 iterations by Realizable two-layer 
k - e (O). Standard k - cu (□), Menter SST k - £J (A) and Spallart-Allmaras (0) models. 



4-1- Angle of attack a — —0.4° 





c x 


c z 


c x error 


c z error 


CPU time, sec 


Experimental results 


0.0419 





n/a 


n/a 


n/a 


Realisable Two-layer k — e 


0.039 


-0.008 


6.9% 


n/a 


225157 


model 












Standard k — uj model 


0.037 


-0.008 


11.7% 


n/a 


203614 


Menter SST k — u model 


0.035 


-0.008 


10.0% 


n/a 


210695 


Standard Spallart- 


0.039 


-0.008 


6.9% 


n/a 


190664 


Allmaras model 













Table 1: The result comparison table for the angle of attack a = —0.4° 
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As can be seen from Table [T] the results for the c z are the same across all 4 
models. The small deviations from the zero value can be easily explained by the 
errors caused by CAD importing and mesh generating procedures, which result 
in small asymmetry of the body. So, in the case of small angle of attack only 
the coefficient c x is compared. As shown in Figures [3]|4] two models, namely 
k — e and Spallart-Allmaras turbulence models, provide the best results with 
the accuracy of 6.9%. The other two models, which are actually standard and 
modified by Menter k — u turbulence models, show much worse results with the 
accuracy of 11.7% and 16.5% respectively. Nevertheless both k — e and Spallart- 
Allmaras turbulence models give very good results, the latter model required 
about 15% less CPU time to finish 1500 iterations (see Figure [2]) ■ It is can 
be easily explained by the fact that Spallart-Allmaras model is a one-equation 
model rather than two-equation k — e model. 

It is suggested that in steady-state simulations of the airships flying with small 
(near-zero) angles of attack, the best turbulence model in terms of accuracy and 
required time and resources is Spallart-Allmaras model. 




i i i i i i i i I a. 

-5 5 10 15 20 25 30 35 <10 ' 

Figure 3: Comparison between the results for the coefficient c x given by the experiment (x) 
and by Realizable two-layer k — e (0)> Standard k — ui (□), Menter SST k — u) (A) and 
Spallart-Allmaras (0) turbulence models. 



4.2. Angle of attack a = 11.62° 

The case of the medium angle of attack gives a bit different picture. In Section 
14.11 it was able to distinguish two best turbulence models in terms of accuracy 
based on comparing the coefficient c x , in this case almost identical results (see 
Tabled for the coefficient c x are presented. 
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Cz 


c x error 


c z error 


CPU time, sec 


Experimental results 


0.0325 


0.233 


n/a 


n/a 


n/a 


Realisable Two-layer k — e 


0.020 


0.248 


38.5% 


6.4% 


231156 


model 












Standard k — w model 


0.019 


0.249 


41.5% 


6.9% 


203065 


Mcntcr SST k — w model 


0.019 


0.249 


41.5% 


6.9% 


228142 


Standard Spallart- 


0.022 


0.247 


32.3% 


6.0% 


190851 


Allmaras model 













Table 2: The result comparison table for the angle of attack a = 11.62° 



The results for the coefficient c z are expectcdly less accurate. This can be easily 
explained by the fact that when the biggest contribution to the coefficient c z 
is made by pressure forces (about 98%), the coefficient c x is mostly (about 
90%) made of shear forces. Pressure forces are much easier to predict than 
shear ones. Thus it is much more complicated to get the same level of accuracy 
for the coefficients significantly impacted by the viscous effects as for the ones 
containing mostly pressure effects. 

So, the result for medium angles of attack gives a very good opportunity to 
compare the accuracy of near-wall layer treatment by 4 different turbulence 
models with the same mesh conditions. Table [2] shows that the best results for 
the viscous effects are shown by the Spallart-Allmaras model. It is remarkable 
that the same model provides the best results for the pressure forces, or, in 
the other words, the highest accuracy level for the coefficient c z . The result 
comparison can be seen in the Figures [3] and 2] 

Taking into account the fact, that the Spallart-Allmaras model uses the least 
CPU time (about 17% less than k — e model, which is the next in accuracy), 
it is evident that Spallart-Allmaras turbulence model is the best model for the 
steady-state airship simulations with medium angles of attack. 



4.3. Angle of attack a = 35.62° 







c z 


c x error 


c z error 


CPU time 


Experimental results 


0.038 


0.825 


n/a 


n/a 


n/a 


Realisable Two-layer k — e 


0.018 


0.895 


52.6% 


8.5% 


236370 


model 












Standard k — w model 


0.011 


0.873 


71.1% 


5.8% 


212138 


Mcntcr SST k — w model 


0.011 


0.907 


71.1% 


9.9% 


212049 


Standard Spallart- 


0.013 


0.890 


65.8% 


7.9% 


183672 


Allmaras model 













Table 3: The result comparison table for the angle of attack a = 35.62° 



For the case of the large angles of attack it is clear (from Table [3]) that with the 
given mesh (especially, in the near- wall layer) all 4 models show quite poor accu- 
racy for the coefficient c x , which is contributed by viscous forces for 75%. Thus, 
the performance in terms of viscous effects prediction is quite poor. In order 



10 



Figure 4: Comparison between the results for the coefficient c z given by the experiment (x) 
and by Realizable two-layer k — e (C0> Standard k — w (□), Menter SST k — uj (A) and 
Spallart-Allmaras (0) turbulence models. 

to increase the accuracy of predicting viscous forces using the given turbulence 
models the mesh in the near-wall region must be refined. 

But the pressure forces have much more significant effect on the overall aerody- 
namic behaviour of the airship, i.e. the coefficient c z have got 98% contribution 
from the pressure forces. Moreover, as the accuracy of predicting the viscous 
forces is almost the same, we can compare the turbulence models based on the 
pressure forces prediction, or, in other words based on the accuracy of the results 
for the coefficient c z . 

Looking at Table[3]one can conclude that the best performance (5.8% accuracy) 
was shown by the standard k — uj model, ft is interesting that the modified 
k — uj models has the worst accuracy of 9.9%. k — £ model requires the biggest 
amount of resources to make the benchmark of 1500 iterations and gives not 
very good accuracy of 8.5%. The Spallart-Allmaras model, which shows the 
best performance in the other two cases (see Sections 14.11 and 14. 2p , again shows 
nice results in terms of resource consumption, but the accuracy is not so good 
in comparison with the standard k — uj model. The comparison between the 
experimental data and the numerical solutions can be seen in Figures [3] and [4l 

So, in case of the large angles of attack, k — uj turbulence model is the best in 
terms of the accuracy, but although Spallart-Allmaras model shows less accurate 
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results it uses about 13% less CPU time for it (for the comparison, see Figure^. 
In this case we can not recommend only one model. The choice of the turbulence 
model should be based on the requirements for the accuracy of the results and 
costs of the calculations. In other words, the choice of the model should depend 
on whether one is ready to trade some extra resources for about 2% of extra 
accuracy, which does not necessarily play a significant role in preliminary airship 
simulation. 



5. Conclusions 



The results presented in Tables [TH3] allow to make a conclusion that among the 
considered turbulence models Spallart-Allmaras turbulence model is the most 
optimal in terms of accuracy (see Figures [3] and 0]) and resource consumption 
(see Figure for the simulations of airships flying at small (near zero) and 
medium (about 10°) angles of attack. 

In case of large angles of attack the standard k — u model performs more ac- 
curate than Spallart-Allmaras model (again, see Figures [3] and 01 but it uses 
significantly more CPU time (see Figure 0). So, Spallart-Allmaras turbulence 
model is suggested as the main choice in the simulation of airships flying at 
large (about 35°) angles of attack as well. 
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Appendices 



Appendix A. Notation for the Spallart-Allmaras model transport equa- 
tion equation 



The turbulent viscosity is expressed as 

v t = vfvi, 

where f v \ = x 3 /(x 3 + an d X '■— vjv. The transport equation is 



(A.l) 
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where 



^ = S -\ 2MJv2, fv2 = 1 — — 7 : 

n l d l 1 + xjvi 
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-C t2 ^[d 2 + g 2 t d 2 t ] 



fw = g 



1 + Cu>3 



<j 6 + ct, 



1/6 



, g = r + C w2 {r - r), r 



Sn 2 d 2 

and d is a distance to the closest surface. 

The constants presented in the equation (|A.2[) have the following values 
K = 0.41; a = 2/3; C bl = 0.1355; C b2 = 0.622; C wl = C bl /n 2 + (1 + C b2 )/a; 

Cti = 1; Ct2 = 2; Ctz = 1.1; Ct4 = 2; C„i = 7.1; C W 2 = 0.3; C w z = 2 
Further details about the model, notation or the constants can be found in Q. 
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